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Abstract 

Using a method of eigenfunction expansion, a stochastic equation is developed for the gener¬ 
alized Schrodinger equation with random fluctuations. The wave field f) is expanded in terms of 
eigenfunctions: b = Yln a n{t)fn{x)i with cj) n being the eigenfunction that satisfies the eigenvalue 
equation iJoAi = where Hq is the reference “Hamiltonian” conventionally called “unper¬ 

turbed” Hamiltonian. The Langevin equation is derived for the expansion coefficient a n (t), and 
it is converted to the Fokker-Planck (FP) equation for a set {a n } under the assumption of the 
Gaussian white noise for the fluctuation. This procedure is carried out by a functional integral, 
in which the functional Jacobian plays a crucial role for determining the form of the FP equation. 
The analyses are given for the FP equation by adopting several approximate schemes. 
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I. INTRODUCTION 


The phenomenon of random fluctuations has been one of the central subjects in statistical 
physics, with topics coveringa wide class of physical systems ranging from molecular level to 
cosmological phenomena [l-4]. Brownian motion, which is one such system, is formulated by 
the Langevin equation, which is an equation of motion that is modified by adding random 
force. The Fokker-Planck (FP) equation, which is written for the statistical distribution 
function for random variables, is an alternative form to describe random systems 1-4], 

The purpose of this article is to study the stochastic analyses of the generalized 
Schrodinger equation in the presence of random fluctuations. The work was inspired by 
previous studies of the Landau-Ginzburg (LG) equation, including the effect of thermal agi¬ 
tation occurring in the superconductivity near the transition temperature (see, for example, 
5j, 6|). We consider this specific problem in the framework of the wider discipline in such 
a way that it can be extended to the general class of the Schrodinger-type wave equation. 
Indeed, among general wave equations, we mention wave propagation in the media in the 
presence of various sorts of fluctuations caused by external agents (see e.g. 7, 8]). Our start¬ 
ing equation is a Schrodinger-type equation with additional random force, which is regarded 
as a Langevin-type equation. Apart from the random fluctuation driven by the external 


agents, the generalized Schrodinger-type equation has been studied in the framework of the 




random potential (9J, [10]. The study presented in Refs. ( 9 , LlOj deals with the randomness 


that is intrinsic to the media itself, and it belongs to a different category from the one that 
will be treated in what follows. 

Although it is easy to write the Langevin equation formally for the wave held -0(a :,t) 
itself, the concrete calculation is rather cumbersome, when one considers applications to 
actual problems. Thus, it is desirable to find a tractable way to apply the calculation. We 
propose here a method of eigenfunction expansion that is based on the use of an orthonormal 
set of functions that are the eigenfunctions of a reference Hamiltonian, Ho, which is taken as 
an unperturbed part of the Hamiltonian appearing in the Schrodinger-type equation. The 
method was previously developed in connection with the semi-classical analyses of a path 
integral jlT , 121. so its application to our stochastic problem is very natural. Therefore, the 
wave held is expanded by this set of eigenfunctions, and the Langevin-type equation for 
the original wave held yields the Langevin-type equation for the expansion coefficients (see 
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the next section). This Langevin equation is converted to the FP equation by using the 
functional integral formalism on the basis of Gaussian white noise for the fluctuation. The 
essential point is that a functional Jacobian appearing in the functional integral plays a key 
role in deriving the FP equation. 

Here, we remark on previous studies of the stochastic approach to the Schrodinger equa¬ 
tion. In Refs. [l3, 14], the stochastic differential equation is employed to describe the 


various sorts of random fluctuations inherent in quantum waves. In particular, in Ref. 


12 |, 


this equation was used to derive the FP equation. The method developed therein may not 
always provide relevant tools for the wide class of problems described by the generalized 
Schrodinger equation. However, the present approach, in which the eigenfunction expan¬ 
sion is used, is expected to provide a very simple and general technique for studying the 
stochastic characteristics of the generalized Schrodinger equation. 

The paper is organized as follows: In the next section, the Langevin equation is derived 
for the generalized Schrodinger equation. In Sec. IIIII the Langevin equation is converted to 
the functional integral on the basis of Gaussian white noise. In Sec. HVi the FP equation is 
derived from the functional integral by adopting the trick of imaginary time. This is based 
on the well-known procedure for deriving the Schrodinger equation from a path integral. In 
Sec. El the analyses of the FP equation are presented for the case in which the expansion 
coefficients are treated as independent of each other. In Sec. IV Al we adopt the strong¬ 
coupling approximation to derive the expression for the “first-excited” state by using the 
variational method. In Sec. IV Bl we present the transition amplitude for the small diffusion 
limit, which corresponds to the semi-classical approximation. 


II. LANGEVIN EQUATION FOR THE SCHRODINGER-TYPE EQUATION 

Our starting point is the Langevin type generalized Schrodinger equation for the wave 
field ^(x,t) jfjjfj: 

7 ^^’ *) ' = t) + V{ X , t), ( 1 ) 

where H is an operator corresponding to “Hamiltonian”. The coefficient 7 denotes the 
complex diffusion constant, and ij(x,t) represents fluctuation coming from thermal and 
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other effects, which obeys the law of white noise, 


(■ rj*(x , t)r](x', t')) = h5(x — x')8(t — t'), 


( 2 ) 


where h is a diffusion constant for the random process. 

In the following discussion, H is assumed to be given by the sum: H = Ho + V = 
(V 2 + f/) + V. H 0 is conventionally called an unperturbative term described by the potential 
U , whereas V represents a perturbation, which allows time dependence in general. H 0 and 
V are assumed to be Hermitian. 

Now let us consider the eigenfunction expansions for 0 and r/, 


i/i(x,t) = y ^a n (t)(f) n (x) , rj(x,t) = J2v n (t)Mx), 

n n 


( 3 ) 


where a n (t ) and rj n (t) are complex functions of t , and 0 n (x) is defined as an eigenfunction 
H 0 (j) n (x ) = A n (ft n (x), where we assume the eigenvalue \ n is positive for the convenience of 
later arguments. By using the orthogonality relation J (p n (x)(p* n (x)dx = 8 nm , we can get a 
Langevin equation for a n (t), 


7- 


da n (t) dE 


+ Vn(t ), 


( 4 ) 


dt da, 

together with the complex conjugate. Here we introduce the energy function E, which 
consists of the unperturbed and the perturbation terms: E — E 0 + E v , with 


E 0 ^ ^ A n a n a n , 

n 

E v = ^2(n\V\m)a* n a m , (5) 

n,m 

where V nm = Vff n (that is, V denotes the Hermitian), and | n) and | m) are the eigenstates 
for mode n and m, respectively. The complex variable a n can be expressed in terms of real 
and imaginary parts as a n = X n + iY n , which is used for the concrete description of the FP 
equation. 

Here we remark on the complex coefficient y in (Hj) : For this purpose, we rewrite Eq. flU) 
in terms of real variables (. X , Y): 


dX r , 


a- 


a- 


dt 

dY„ 

dt 


+ K- 


dY n 


dt 

dXr, 


k- 


dt 


1 dE 

2 dX^ 
1 dE 

~2dYr 


+ Re(r/), 
+ Im(?/). 


( 6 ) 
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Here a and k are defined by 7 = —a + in. If we set k — 0, this equation describes the 
dissipation, hence a should be positive. On the other hand, for the opposite case in which 
a = 0 , one can see that the equation of motion leads to the non-dissipative motion, since 
the energy function E is conserved. 


III. FUNCTIONAL INTEGRAL 


According to the assumption of white noise, rj(x,t) obeys the Gaussian distribution; 

1 

2h 


P [rj(x , t)] = iVexp 


rf (x, t)r](x , t)dxdt 


(7) 


where h is the diffusion constant corresponding to the Planck constant. By considering 
the Langevin equation and taking account of the boundary condition: -0(x,O) = ^(x), 
?/>( x , T) = ^'( x ), the transition probability from x ) to is written by the path integral, 


/■b'O) 

K(^'(x), T\ijj(x) } 0) ~ / exp 

J \jj(X ) 


2 h 


rj*(x, t)rj(x, t)dxdt 


II Vr){x,t)*Vr)(x,t), ( 8 ) 

X,t 


which turns out to be represented in terms of the amplitude {a n }, 


r{<} 

K ({ a 'n}, T \{ a n},0)= exp 

^ {^n} 


■— r 

2ft Jo 


^2Vn(t)Vn(t)dt 


Yl Vr ln(tyVr]n(t), (9) 


n,t 


where we use the eigenfunction expansion defined by Eq.(151h The functional integral 


is 


rewritten in the following steps: First, by inserting the expression for the identity of the 
(5-functional integral |15|, [l 6 j: 


[ Unit)- 7 

't n t ' 


da„ dE 


( 10 ) 


d t da* 

together with use of the well-known relation for <5 function: 5[f(x)] oc Jf^exp [ipf(x)\ dp , 
and next by integrating over r/ n ( 77 * ) as well as p n ( p* ), we obtain 

r- T 


K ~ C / exp 


'2 hj 0 T, F n Fndt 


VF:VF n , 


( 11 ) 


where C is a factor that does not depend on F n . By the definition of F n and by using a 
functional determinant, the path integral is given as 


K ~ C / exp 


1 

2h 


j2(i* da * n 


dE 


dt da r , 
xdet 


da n dE 

1 


Va* n 'Da n . 


r<5F*i 




det 


U«nJ 

_ 3cin _ 


( 12 ) 
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The functional determinant is calculated and written by using real variables (A", Y) 


det 


\dFY 


r 1 r T 

f d 2 E 

Q2E \ j 1 

_ n _ 

= exp 

. 87 y 0 

Kdx 2 n 

+ by;-) dt \ 


( 13 ) 


This factor plays a crucial role in fixing the correct form of the FP equation. The process 
to derive it will be given in the Appendix [A] 

From the above form of path integral, we see that an action functional is modified as 


s = E 

n 

= E 


T 


da ' as V|(l + i)E 


dt dat 


" T / d 2 E t d 2 E\ u 

ax! + W 2 dt 


r o 


£ n dt + h 


Mdt 


(14) 


where 


M = 


a 


d 2 E d 2 E 


(15) 


4|tI 2 \8Xl 3Y,l 

Now we rewrite the Lagrangian C n in Eq. (ITT) by using real coordinates, (A, Y), explicitly, 
and further by noting 7 = —a + in. Then we obtain 

1 


= ^<h\ Z (X 2 + Y 2 )+a(X ngx 


dE ■ dE 
+ Y n 


-nix:M-Y„ dE 


dY n 


dX r 


+ 


dE 

dXl 


dY„ 


+ 


dE 

W n 


(16) 


Here in order to make a connection with quantum mechanics, we introduce imaginary 
time r(= —it), and we set the Lagrangian as C n {— — £ n ). By using the conjugate momentum 


Pn = Iw 3 - (X n = l j jl ), the Hamiltonian corresponding to C n becomes 


U n = 


2|tI 2 

1 


PnX ~ IOL 


dE . dE 


dX r 


+ IK 


dY n 


dE dE 

+ PnY ~ - IK- 


dY n 


dX r , 


dE 

dXr 


+ 


dE 

w n 


(17) 


Thus, the total Hamiltonian is given by 'H — %, 


IV. THE FOKKER-PLANCK EQUATION 


We now derive the FP equation via the “imaginary time” Schrodinger equation by starting 
with a quantized version of the Hamiltonian (ITT)) . First we note that the corresponding 
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“wave function” \l/({X n },T) is related with the probability distribution function P({X„},f) 
as follows: 


*({X n },r) = P({X n },t). (18) 

As is well known, the Schrodinger equation is given by integral equation 
^({X„},r + e)= f K({X n },r + e | {XJJ, r)V>({X'„}, r) 

J n 

where the propagator K is constructed from the action functional given above. Thus, accord¬ 
ing to the well-known procedure (see, 

T({X n },r) turns out to be 


17]), the Schrodinger equation for the wave function 


^ 1 

lh ~dr = ^2hf 

n 11 


a dE k dE \ f a dE k dE 

VnX-1-^TT+l- 77 ^ + 


2 dX n 


2 dY n 


2 dY h 2 dX„ 




(19) 


where 


PnX = -ih 


d 


W = 


dX, 
dE 


PnY = ~ih 


()X„ 


d 

dYn 


dE 

w n 


+ Mh. 


( 20 ) 


Then, by using Eq. ([20j) . and replacing the imaginary time r with the original real time t, 
namely, t —it, we obtain 


2|t|- 


a 


dE <9T 


<9T 

~dt 

+ 


E- 


h 


n 

dE <9T 


dX n dX n dYn dY r 


2| 7 | 2 

+ K 


d 


+ 


d 


dX n ) \dY r 

dE <9T dE 




dX n dY n dY n dX r , 


( 21 ) 


a 

(d 2 E 

d 2 E\ „ J 
+ M 

UiV 


W 2 J J 




Using ffl5]h the last term becomes ^yp Hence the second and the last terms 

can be combined in a single form by noting the relation VA-/ + A-V/ = V - (A/). Further, 
by getting the wave function T back to the original probability distribution P, we finally 
obtain the FP equation: 


dP 

dt 


E 


h 


2| 7 | 2 \dX , 


d 


h 


d 

Wn 


P + 


P + 


1 dfdE dE „ 
a——P — k^t^P 


2| 7 | 2 dX, 

1 d 


2H 2 dY n 


dX n 
dE 


dY n 


a 


dYn 


dE 

P + K 3X/ 


( 22 ) 


















































Here we discuss some specific features of this equation: 

(i): Equation (l22j) can be represented in a continuity equation of the probability: 


dP 

a +v ' J = °- 


where the “current” is defined as J = Yh n Jn, 

1 


J n. 


,. _ „ . . d E d E 

( hV n P + A n P) , A n = a^z -kkx 


( 23 ) 


(24) 


2| 7 | 2V ""‘ “ax„ .ax„’ 

with k being the unit vector that is perpendicular to the (A, Y) plane. Indeed, it has been 
known that the FP equation can be derived by following an analogy with the continuity 


equation for the probability flow in an intuitive way 
(ii): As a special < 
for which we propose 


18] 


ii): As a special case, we consider stationary distribution, namely the case of ^ = 0, 


P = e- 0E , 


(25) 


where (3 = 1/kT denotes the inverse temperature. Substituting this form, we have 


(a 


- w) E 


d 2 E d 2 E 

+ 7W7TT-^ 


dK 


dY 2 

n 



= o. 


(26) 


For this relation to be satisfied for arbitrary function E, the following relation should be 
held: 


P = 


a 

h' 


(27) 


This is merely the fluctuation dissipation theorem, which establishes the relation among 
three parameters a, /3, and h, representing the dissipation, inverse temperature, and dif¬ 
fusion, respectively. Note that the effects of n disappear in the above relation ([26]) . and 
non-zero a plays an essential role for the existence of the equilibrium state. If a = 0, the 
equilibrium state does not exist, because the effective temperature becomes infinite. 

Calculation of the partition function : The equilibrium distribution is used to evaluate the 
partition function, which leads to various sorts of thermodynamical quantities [ll|: As a 
particular case, we consider the case in which the energy function is given by the unperturbed 
one, E 0 = \ n (X 2 + Yf), for which we have 


Z 


exp[— j3 A„a*a„] da* n da n 


n 


2vr 

Wn 


(28) 
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This is just the same form as the partition function used in the superconductivity near 


the phase transition jf], 20], which is described by the time-independent LG free energy. 
That is, the partition function is given by Z = f exp [— (3F [ip, ip*]\ ]Q dip*dtp with F[ip, ip*] = 
f ip*Hipdx. H is given as the Hamiltonian that depends on the electromagnetic field. In this 
way, the present approach includes the stationary problem for the LG theory as a special 
case. 

(iii): The FP equation describes the coupling among the modes a* (i — 1, ■ • ■ , oo). This 
is governed by the perturbation term E v in the energy function (l5j) . which is written in terms 
of the real variables (X n , Y n ): 


E v — B nm (X n X m + Y n Y m ) + C nm (X n Y m — X m Y n ), 


(29) 


n,m 


where H, C are the real and imaginary parts of the hcrmitian matrix V nm : ( V nm = B nm + 
iC nrn ); hence these become the symmetric and skew-symmetric matrices respectively, namely, 

Bnm B mn , C nm C mn . 

In particular, the FP equation described by the unperturbed energy Eq yields 


_ v 


h 


2| 7 | 2 \dX, 


d 


h 


d 


P + 


P + 


An d 
\l\ 2 dX, 

An d 


(aX n P - kY u P) 


(aY n P + kX v P) 


(30) 


2| 7 | 2 \dY n J | 7 | 2 dY n 

Furthermore, the term coming from E v gives rise to the modification to the above; that is, 

1 d 


E 

nm 

E 


h\ 2 dx r 

1 d 


1 2 dY n 


« (BlmXm + Ci m Y m ) P — K ( Bi m Y m — Ci m X m ) P 
& (BlmYn Cl m X n j) P T K ( Bi rn X m ClmYmP) P 


(31) 


This term plays the role of the perturbation to the unperturbed equation (1301) . Also the 
coefficients B nm , C nm have time-dependence in general, so we have a time-dependent per¬ 
turbation theory for the FP equation that may be expected to bring about specific physical 
consequences. 


V. ANALYSES OF THE STOCHASTIC EQUATION 

The FP equation that was derived above is described by the energy function E : which 
includes the coupling among various modes. In the following discussion, we restrict the 
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argument to the specific case of the independent mode E 0 = For this case, 

the stochastic equation is given by (l30]h and the analyses of this will be carried out in a 
simple manner, because the modes are decoupled from each other. That is, the distribution 
function is written as a product form: 

P({X} y t) = l[P n (X n ,Y n ,t), 

n 

where P n (X ni Y n ,t) means the distribution for the n-th mode, and the term including k is 
omitted because it does not contribute to the equilibrium state. In what follows, we consider 
the n-th mode only, so the index n is suppressed. Then, we set the variable separation as 


P(X,Y,t) — exp[—et]f(X,Y). 


Thus, we get the eigenvalue problem as follows: 


/ d 2 


+ 


d 2 

dY 2 


f + 


2A a 
~h~ 


x^- + 

dx dY 


‘In 

~(m + 2X)f. 


(32) 


Here we use the polar coordinates: A" = R cos 6, Y = R sin#, and we consider the case in 
which the eigenfunction f(X,Y) does not depend on the angular variable 6. Therefore, the 
eigenvalue problem becomes 


IA \ B (sL + aMf) 

RdR [ \dR + 2 dR J ) 



(33) 


where ji = = 2A/3. Here we use the notation U(R) = X 2 + Y' 2 = R 2 , then Hamiltonian 

is written as H = XU(R). 


A. Strong coupling approximation 

The eigenvalue problem (l33]i looks simple, but this may not be represented by special 
functions that have been used so far. If we note that Eq. (j33T) includes the parameter //, 
which can be regarded as a perturbation parameter, one may think of carrying out the 
perturbation scheme by expansion with respect to the small parameter /i (|p| 1). That 

is, one starts with the equation for p = 0 as an unperturbed solution, which is given by 
the Bessel function. However this procedure may not be relevant, since the case (i = 0 
corresponds to a = 0, which does not represent thermodynamic equilibrium as it is pointed 
out above. From this inspection, it is suitable to consider the case in which \fi\ is not small, 
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~ 1, which we call the “strong coupling approximation”. The following procedure is 


similar to that used in the stochastic approach of ferromagnetic particles 
problem under consideration belongs to a completely different discipline. 
Then, we immediately get the solution for the zero eigenvalue (e = 0): 

MR) = [~U(R) 


18j . although the 


(34) 


This corresponds to the Boltzmann distribution, which is identified with the “ground state”. 
As a sequel of the ground state we assume “excited state” as 

ho 


f(R)=ex P ~U(R) g(R). 


The function g(R) satisfies 


d 

dR V 


( Rex p 


h 


U 


dg 

dR 


2ch 


ei?exp 


h 


U 


(35) 


(36) 


In particular, for e = 0, we can set ^ to be zero, which gives go(R) = k (constant). 


The eigenvalue equation 
the functional given by 


is equivalent to the following variational problem: That is, 


/ = 


'0 


dg 

dR 


exp 


h 


U 


RdR 


should be minimized under the normalization condition 

/•no 

h, 


g 2 exp 


'0 


-U 


RdR = 1. 


(37) 


(38) 


Furthermore, we propose that the orthogonality relation should be held between two eigen¬ 
states g and g\ 


gg exp 


-»U 


RdR = 0. 


(39) 


Under this prescription, we consider the first-excited state as a concrete example. To 
perform the variational calculation, we choose the trial function as the quadratic expression 


gi (R) = A + BR + CR 2 . 


(40) 


The coefficients A, B and C are determined by three conditions corresponding to Eqs. (1371) - 
(j39|) . which are written explicitly as 

I = RB 2 + 4 J 2 BC + ARC 2 (= e,) (41) 
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and 


RA 2 + 2 J 2 AB + J 3 (2AC + B 2 ) + 2 RBC + J 5 C 2 = 1, 

JiA + JoB + t/sC* = 0. 


( 42 ) 


respectively. Here we set Jk = / 0 °°exp [— R k dR. The minimum value of / is obtained 
as follows: First, by eliminating A from Eq. (l42l) . one gets the quadratic constraint with 
respect to B, C\ 

G(B , C) = + J 3 ) 5 2 + 2 RBC + (j 5 - C 2 - 1 = 0. (43) 

Next, by using this constraint and following the Lagrange multiplier method, we have the 
relations 

JL(I-\G) = 0 ( Xi = B,C, A), (44) 

where A is a multiplier. Equation iITTli leads to simultaneous equations for B , C, and A. By 
solving this, we can obtain the value of e\. The explicit manipulation is rather tedious and 
is omitted here. 

By continuing the above procedure, we can obtain a sequence of the excited states g n (R)- 
By using this, the distribution function P(R , t) is expanded as 


P(R,t) = C n (t)g n (R) exp 




(45) 


The expansion coefficient C n [t ) satisfies dC n /dt = —e n C n , and from the initial condition 
P(R, o) = Y^n C n(0)g n (R) exp [-f U], we get C„(0) = / P(R, 0)g n (R) exp [-f U] RdR. As 
an actual situation, it may be sufficient to keep the terms up to the first excited state, 


P(R,t) = {C 0 («) 9 o + Ci(*) 9 i(^)}exp 



(46) 


The distribution of this form may be utilized to calculate the time evolution of the average 
value of physical quantities under consideration. 


B. Small diffusion limit 

We consider another procedure of the nonperturbational scheme that is given as the 
asymptotic limit of the path integral, namely the limit of zero diffusion, h ~ 0. This 
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corresponds to the semiclassical approximation in quantum mechanics. The path integral 
expression for h ~ 0 becomes 


K 


SC 


exp 



(47) 


Here the classical action is written as a sum of the contributions coming from each mode rt; 
S sc = Xn Ssc- I n what follows, we take account of the mode n only, and we suppress the 
index n. The hnal result is obtained by summing over n. By using the polar coordinate, the 
action functional is given as 


S„,= 


a 

T 




dt, 


(48) 


where ( = The conjugate momentum for 9, pe = R 2 9 , is conserved because 9 is a cyclic 
coordinate. Then, the Lagrangian becomes 


L = 


a 


R-C 


dU 

dR 


+ 


2R 2 


(49) 


Following the procedure of classical dynamics (see [2l|), the equation of motion is derived 


by a Routhian given by 1Z = c9 — L. ft is rather complicated to deal with this equation, 
but one can ignore the term including c if c can be chosen as small. Then the equation of 
motion is written in a form of the so called “instanton” equation 

dU 




(50) 


By using this, K sc is simply written as 

K sc = exp 



(51) 


Furthermore, noting = 7 ^, it follows that the transition probability from the initial 
point Ri to the final one Rf is 


K sc (Rf,t\Ri,ti) = exp 


r c 2 

r R f dR 

1 

JRi ^ 3 J 


= exp 


8 Ch \ R\ Rf 


(52) 
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As the final step, attaching the index of mode for the orbit, we write R as R(k), and 
taking the product over all modes k, we obtain 


K sc ({Rf},t\{Ri},ti) = JJexp 

n 


4 

[ R ( k) f dR(k) 

4C kh 

!R(k)i R(k) 3 


= U exp 
k 


8 (kh \R{k )? R(k)} 

where c*,, Cfc indicate that the quantities correspond to mode k. 


(53) 


VI. SUMMARY 

We studied a stochastic theory for the generalized Schrodinger equation by using a method 
of eigenfunction expansion. The present approach would have an advantage in that once one 
starts with the expansion for the wave function (order parameter) in terms of the set {a n }, 
one can always obtain the Langevin equation in a very general way. Thus this approach 
would provide an efficient way to investigate a wide class of systems that can be described 
by the generalized Schrodinger equation. 

By using the functional (or path) integral formalism, the Langevin equation results in 
the FP equation based on the assumption of Gaussian white noise for the fluctuation. We 
paid particular attention to the calculation of the functional Jacobian, which is simply 
incorporated in the action functional. As a consequence of this procedure, we arrived at the 
specific form of the FP equation (l22lh which is relevant to our subsequent analyses of the 
distribution function. 

The analyses of the distribution function have been given for the case of the unperturbed 
Hamiltonian, which can be treated within two categories: (i) by using the expansion with 
respect to the dissipation constant a, and (ii) by adopting the asymptotic limit of “zero” 
diffusion {h ~ 0). Our analyses have an advantage in that we can obtain a concise analytic 
form for the distribution function and the transition probability. In this connection, the FP 
equation that was developed in [p| seems to be a complicated way to obtain a simple form 
of the distribution function. As a potentially useful application of our method, we mention 
the problem of a calculation for various sorts of transport coefficients in nonequilibrium 
statistical physics that is inspired by the generalized Schrodinger equation 22]. 

An immediate application of the present formalism could be the calculation of conduc- 
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tivity for the superconducting fluctuation current just above the critical temperature. The 
FP equation (|22jl is expected to play a key role in calculating the time-dependent average 
for the current leading to the conductivity. The details of this topic will be discussed in a 
forthcoming paper. 


Appendix A: The calculation of the functional Jacobian 


To evaluate det 


6F n 

Sdn 


, we start with a general procedure that was developed in 23]. Let 


us write F = 7 ^ — A, and define F = 7 (A We note the relation 


db 


dt 


dt' 


det 


FF 

5a 


= det 


5b5b 
5b 5a 


= det 


5b 

5b 


det 


5b 

5a 


(Al) 


The first factor det 


= det [A] — Q can be omitted because it is simply a divergent 


factor. Then we have the integral equation 

lfb{t) = 70 (f) - 


A [a(r)] dr, 


(A 2 ) 


which is rewritten by using the step function, 

76 (f) = 70 (f) — f 9(t — t)A [a(r)] dr. (A3) 

Jo 

Using this form, the interval of integration can be converted to [0,T], which coincides with 
the interval in the path integral (fT2jh Then, Eq. (1 A> I) is expressed as a discrete form: 


7 b{ti) = 7 o(f;) - ^ Ofa ~ T i) A KU')] £ - 
3 

By carrying out the differential with respect to a(ffc), we obtain 

db(tj) _ da(tj) _ _ dA [o(rj)] 

7 da(t k ) 7 da(t k ) ^ 1 J da(t k ) 

which is reduced to 


(A4) 


(AS) 


dbiti 


= 5i 


da(t k ) “ lk 7 " da(t k ) 

The matrix represented by these elements is a triangular matrix. For t t > t k , it follows that 


-0(U 

7 


t K ) 9A [o(f*)] c 


(A 6 ) 


db(tj) _ 1 dA [g(tj )] 

da ( 4 ) 7 da(ti) 


(A7) 
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For ti = tk, noting 9{ti — tk ) = 1/2, which is a Dirichlet discontinuous factor 24], we get 


db(ti) 


= 1 


1 dA[a(U)\ 


e. 


(AS) 


da(U) 27 da(U) 

Here it is crucial to have the factor 1. The determinant, except for the infinite factor C, is 
given as 


det 


5_F_ 

5a 


= det 


5b 

5a 


n> 


1 &4[a(t,)] 

27 da(U) 


which turns out to be 


det 


5b 

5a 


= exp 


L 27 


dA[a(t)] 
da(t ) 


dt 


Therefore, by substituting A = J^r, we have the functional Jacobian, 


det 


rj/Ai 


r 1 r T 

fd 1 2 3 4 5 6 E 

d 2 E\ 1 

_ Sd n 

= exp 

1 

1 

001 

-51 

0 * 

[dx* 

+ w —0 I dt 

dY 2 ) \ 


Also we have the complex conjugate of (IA 11 D . 


det 


5F* 

w n 

5a* 


= exp 


/ d 2 E 8 2 E . , 

I —— + ) dt 


[ 87 * Jo \dx 2 oy;? 


(A9) 


(A10) 


(All) 


(A12) 


By summing up (IA1 1|) and (1 A 1 2 | 1 . we obtain the terms that are proportional to M in the 
text. 
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